The Vector Representation of Standard Quantum Process Tomography 



Xiaohua Wu 

Department of Physics, Sichuan University, Chengdu 610064, China. 

The characterization of the evolution of a quantum system is one of the main tasks to accomplish to 
achieve quantum information processing. The standard quantum process tomography (SQPT) has 
the unique property that it can be applied without introducing any additional quantum resources. 
In present work, we shall focus on the following two topics about the SQPT. At first, in the SQPT 
protocol for a d-dimensional system, one should encounter a problem in solving of a set of 
linear equations in order to get the matrix containing the complete information about the unknown 
quantum channel. Until now, the general form of the solution is unknown. And a long existed 
conviction is that the solutions are not unique. Here, we shall develop a self- consistent scheme, in 
CN| , which bounded linear operators are presented by vectors, to construct the set of linear equations. 

With the famous Cramer's rule for the set of linear equations, we are able to give the general form 
c I ^ ■ of the solution and prove that it is unique. In second, the central idea of the SQPT is to prepare 

04 \ a set of linearly independent inputs and measuring the outputs via the quantum state tomography 

^ . (QST). Letting the inputs and the measurements be prepared by two sets of the rank-one positive- 

^ ■ operator- valued measures [POVMs], where each POVM is supposed to be linearly independent and 

I informationally complete (IC), we observe that SQPT now is equivalent to deciding a unknown state 

with a set of product IC-POVM in the d'^-dimensional Hilbert space. Following the general linear 
state tomography theory, we show that the product symmetric IC-POVM should minimize the mean- 
square Hilbert-Schmidt distance between the estimator and the true states. So, an optimal SQPT 
can be realized by preparing both the inputs and the measurements as the symmetric IC-POVM. 
keyword: Standard quantum process tomography; The Cramer's rule; Linear state tomography. 
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qH. I. INTRODUCTION 

The characterization of the evolution of a quantum system is one of the main tasks to accomphsh to achieve quantum 
^ ■ information processing. A general class of methods, which have been developed in quantum information theory to 
lO , accomplish this task is known as quantum process tomography (QPT)- for a review of quantum tomography, see Refs. 

■ [1 — 3] . The various protocols of getting the complete information about the quantum process can be divided into two 
' } \ classes: The standard quantum process tomography (SQPT) [1,4,5] which does not require any additionally quantum 

■ resources, and the so-called ancilla-assisted quantum process tomography [6-8]. 

[ Usually, the complete characterization of x matrix, which contains the complete information about a unknown 
quantum chnnel, is a non-scalable task: For the N c?-level system, there are about d^^ elements to be decided. 

■ Recently, a series of works, where the ancilla system is also required, have demonstrated that it is possible to extract 
. . , partial but nevertheless relevant information about the quantum process in an efficient and scalable way [9-15]. These 
^ ' approaches share an essential feature that: They are based on the idea that the tomography of a quantum map can be 

k>( , efficiently performed by studying certain properties of a twirling of such a map. Another method, the so-called direct 

■ characterization of quantum dynamics (DCQD), was also constructed for use in partial characterization of quantum 
d \ dynamics [16-17]. In experiment, the selective and efficient quantum process tomography has been recently carried 

out with single photons [18]. 

With the aid of a clean ancilla system, the approach of ancilla-assisted quantum process tomography is to encode 
all information about the transformation into a single bipartite system-ancilla quantum state, and thus completely 
reduce the problem to that of quantum state tomography. In general, there three known methods of performing QST: 
The maximum likelihood [19-20], the Bayesian mean method [22-25], and the linear state tomography [26-27]. Our 
present work is motivated by such a simple reasoning: To carry out the ( complete and partial ) ancilla-assisted QPT 
in experiment, one should at first perform a complete QPT experiment to show that the ancilla system is (nearly) 
perfect. To avoid the logical recycling, this complete QPT experiment should be realized by the SQPT protocol which 
does not require any additional resources. In this aspect, the SQPT is not only the unique method, in which the 
quantum process tomography can be directly performed, but it also plays an important role in realizing the various 
kinds of ancilla-assisted SQPT in experiment. 

In present work, we shall focus on the following two topics about the SQPT. At first, in the SQPT protocol for a 
d-dimensional system, one should encounter a problem in solving of a set of d'* linear equations in order to get the 
matrix containing the complete information about the unknown quantum channel. Until now, the general form of 
the solution is unknown. And a long existed conviction is that the solutions are not unique. Here, we shall develop 
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a sclf-consistcnt scheme, in which bounded hnear operators are presented by vectors, to construct the set of hncar 
equations. With the famous Cramer's rule for the set of hnear equations, we are able to give the general form of the 
solution and show that it is unique. In second, the central idea of the SQPT is to prepare a set of linearly independent 
inputs and measuring the outputs via the quantum state tomography (QST) . Letting the inputs and the measurements 
be prepared by two sets of the rank-one positive-operator- valued measures [POVMs] , where each POVM is supposed 
to be linearly independent and informationally complete (IC), we observe that SQPT now is equivalent to deciding 
a unknown state with a set of product IC-POVM in the d^- dimensional Hilbert space. This observation makes it 
possible for us to judge which choices of the inputs and measurements are optimal for SQPT by following the general 
discussion of Scott [26-27]. 

The content of present work can be divided into following parts. At section 2, we shall introduce the isomorphism 
where bounded operators are related to vectors. Several super operators and some relations among them should be 
carefully defined and proven there. In section 3, we shall represent the SQPT theorem in the vector denotation. The 
general solution of the set of linear equations should be given. Its uniqueness should be proven by the well-known 
Cramer's rule. In section 4, by jointing our SQPT protocol with the linear state tomography theorem, we show 
that the symmetric IC-POVM, which has been introduced [28], should be used as the inputs and the measurements 
operators of SQPT. Under the linear quantum state tomography, the biproduct SIC-POVM not only feature a simple 
state reconstruction formula but also minimize mean-square Hilbert-Schmidt distance between the estimator and the 
true states. In section 5, for the optimal case, we offer a simple formula to get the matrix containing the complete 
information of the quantum channel. Finally, we end our work with a short discussion. 



II. THE VECTOR DENOTATION 



Consider a d-dimensional Hilbert space H with {\i)}f^i the basis for it. We could define a set of operators {Cij}f j^i, 

Cij = \i){j\,Tr[ClCkj]=5ikSji. (1) 

This set of operators formulate a orthogonal basis on the state space. Any dxd matrix in H, say A, can be expanded 
with it, 

d 

A = AijCij,Aij = Ti[CljA], (2) 

with Aij the matrix elements of A. In the following argument, we will make \mc of the convenient isomorphism in 
which bounded linear operators are presented by vectors and denoted with the symbol [•)). It should be noted that 
the ways of defining [•)) are not limited [29-33]. Although our scheme is similar with the one introduced in [33], for 
the self-consistency of present work, it will be presented in an independent and complete form. 
Theorem 1. Let \S+)) be the maximally entangled state in iJ®^, 



\s+)) = ^Y.\^)®\k). 



If an isomorphism between the d® d matrix A and a d^-dimensional vector \A)) is defined as 

\A)) = V~dA®la\S+)), (3) 

there should be 

d d 
|A)) = ^ ^ (4) 

Proof. From the definition of the isomorphism above, one may check that each operator in (1) is related to the basis 
vector in iJ®^, 

\Cij)) = ii> ® \3) ^ m- 

Since that the isomorphism is linear, we shall get the form in (4) where the matrix elements A^j are related to the 
expanding coefficients of the vector \A)). 
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Theorem 2. Suppose that A , B, and p are three arbitrary bounded matrices in H, there should be 

Tr[A^B] = {{A\B)), (5) 
\ApB))=A<E>B'^\p)). (6) 

Proof. The equation (5) can be proved with the definition of the vector in (4), 

{{A\B)) = ( J2 ^ijmv E ^fc'i^o) = E4^^^- = 

k.l=l 

The equation in (6) can be verified in the way fike, 

d d d 

A(g)B'^\p)) = {J2 Ak\i){k\)(E>{Y,Bij\j){l\) J2 Pmn\p)) 

i,k+l j,'=l m.n=l 

d 

= E ^ikPklBlj\ij)) 

i,j,k,l — l 

= \ApB)). 

Let the quantum channel of the H be described with a set of Kraus operators, e : {E^}. With p for the input of 
the channel, the output £{p) should be decided by the Kraus operators, 

m m 

The channel is supposed to be trace preserving. For this set of Kraus operators, a supper operator is defined as, 

X^ = Y1 \En){{En ^d-e(^ U\S+)){{S+\), (8) 

m 

which is Hermitian and positive. Its matrix elements, Xij^ki ~ ((^ilx^l^O)' ^^^^ form 

Xlr,ki=T.^^™*- (9) 

m 

If these coefficients have been known, the output in (7) should be decided with the given input, 

<P) = J2xtr,kS)ijm{k\. (10) 

ij;kl 

Noting Tr[x*^] = d, we can also define a normalized density operator pe^ 

p, = ^x',Pe = e®U\S+MS+\). (11) 

Prom it, we see that the tomography of x'^ is equivalent to determining the unknown state p^ in H'^'^. With {\ij))}fj^i 
the basis of if®^, the set of operators {C'ij;/si}f,j,fe,;=n which are defined as 

Cir,ki = \ijMkl\ = \i){k\^\j){l\, (12) 
should formulate an orthogonal operator basis. Any arbitrary cf (g) cf matrix, say, F, can be expanded with it 

d 

r= E ^irM^irM,'^ij;kl =T^[Cjj.kl^], (13) 
with Tij-ki the matrix element of F, Tij-^i = {{ij\T\kl)). Jointing (8-12) together, we shall find the result 
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and rewrite (10) with 



e{p) = Trp)(i| ® \l){j\pM{j\p\m 

ij;kl 



where the so-called Jamiolkowski isomorphism is recovered here [34]. 

As a direct application of (6), we find that the two vectors, \s{p))) and \p)), are related by 

m 

For convenience, we could introduce a super operator A'^, 

A= = E£;™®(£;™)*, (14) 

m 

and let the two vectors, which corresponds to the output and input of the quantum channel respectively, be related 
with a simple formula, 

|£(p)))=A1p)). (15) 

For the same set of Kraus operators {E"^}, we have introduced two different super operator, x'^ in (8) and A'^ in 
(14), here. There is a question: In which way these two super operators are related? To answer this question, we 
shall at fist define another isomorphism where a. cP ^ (P matrix F is related to the vector |r) in H^'^. 
Theorem 3. Let |5+) be a maximally entangled states in iJ®^, 

we could define a vector |r), 

\T) = d-T()^laAS+). (16) 
With Fjj.fe; the matrix elements of F, Fj^.;;; = {{ij\T\kl)), there should be 

d 

|F)= E YirM\i3\kl),\ir,kl) = \ij))®\kl)). (17) 

Proof. It can be shown that each basis operator Cij-ki in (12) should be related to the basis vector kl), 

\Cij;ki) = \ij;kl) = ^ \kl)). 

With the known fact that the isomorphism is linear, we can get (17) by using the result in (13). 

It should be noted that the isomorphism in (16) differs from (4) only in its dimension of the Hilbert apace. So, we 
may introduce the following theorem without giving proof for it. 

Theorem 4- Suppose that F, A, and S are three arbitrary bounded matrices in H^^, there should be 

Tr[FtA] = (F|A), (18) 
|FSA) =F0 A'^IS). (19) 

As the second step, we shall introduce a super operator (3^ in iJ®^ with the following theorem: 
Theorem 5. Let 13'' be an operator in i/®^, 

d 

E \ij;kl){ik;jl\. (20) 

i,j,k,l=l 

It is Hermitian and unitary, 

0^ = (/3^)t = (r)-i. (21) 
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The two vectors, \\A)){{B\) and \ Aigi B*), are related by 

m){{B\) = \A^B*). (22) 
Proof. The property of ^'^ in (21) can be easily verified with the definition for it. Equation (22) can be proved as 

i^j.k.l — l 
d 

d 

= E Ai^BU\i){3\®\k){l\) 

i,j,k,l 

= \A^B*). 

Finally, as a direct conclusion based on the results above, the answer of our question is known: The two super 
operators, and A'^, are simply related by the equations in below, 

|X'=)=^1A^),|A^)=/3^|X^). (23) 

After having defined the two super operators, x'^ and A*^ above, for describing the quantum channel, we shall 
continue to introduce the operators which are used for the measurement. Suppose that {Pf,}f^^i is a set of linearly 
independent operators ( which are not necessary to be orthogonal with each other). According to the isomorphism in 
(4), wc shall get a set of linearly independent vectors {|-P/j))}^=i- Following the idea of Chefles [35,36], it's convenient 
for us to introduce the reciprocal state \Rij,)), which is defined as that lies in iJ®^ and is orthogonal to all \Pv)) for 

= V. (24) 

Such a set of reciprocal states {|-R^))}^Li always exist if and only if the states in {|-P^))}^Li are linearly independent 
[35,36]. Before proceeding with the further discussion about these operators, we shall at first make our denotation 
into a more compact form. In present work, we suppose that all the calculations in iJ®^ are carried out with the basis 
{Ki))}ij=i- ■f'or simplicity, we may introduce the denotation, 

]M)) = |^j)),/x=(rf-l)^ + i, (25) 

where the single Greek letter /i takes values from 1 to d^. Therefore, one may use F^;^ for the matrix elements of F. 
In the calculation, Fj^.y = {{ij\r\kl)), the actual values of the indices, i, j, k, and /, should be decided by the rule in 
(25). Now, let's return back to the two set of states in (24) and begin at the following theorem. 

Theorem 6. For the set of linearly independent operators {P^j^^i, the super operator P is defined by, 

P = Y,\P^}}mdet{P)j^O. (26) 

It has the properties that 

P\l^)) - K)), (27) 

\R,)) = {p-y\n)), (28) 

j2\R^))m='^d^ (29) 

Proof. The first equation can be get with the definition of P. Since that det(P) 7^ 0, it's inverse P^^ is unique 
decided by P~^ = ^^^i Based on it, the equation (28) can be verified. The final result can be also verified 

with (27) and (28): E^i = Ef=i(^-') VX^I^^^ = i^,. 
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In fact, the set of reciprocals states introduced above can be also defined as the canonical dual frame corresponding 
to {|P^}f=i. 

Theorem 7. For a known set of linearly independent vectors {|-P^))}^=i, a super operator is defined by 

:F=J2\Pf^)){{Pf^ld^iF^0- (30) 

The canonical dual frame corresponding to {Pij.}^^i can be defined as 

\R,))=J^-'\P^)). (31) 

From these definitions, there should be ((i?y|P^)) = (5,^^. 

Proof. Using the definition of the super operator P, one may get T = PP^. Since that detF ^ 0, the inverse of 
F exists in an unique way, = {P~^y P~^ . Jointed it with (27), there should be 

:f-'\r,)) = (p-i)tp-ipiM)) = (P-^)V)), 

where the reciprocal states defined in (28) is recovered here. 

For the special case where the set of operators {Z)^}^^i are orthogonal with the algebra Tr[Dj^Di,] = dfj^v, one may 
directly define a super operator U, 

U = Y,\D,))m. (32) 

It can be easily verified that U is unitary since that 

U\H)) = \D,)) (33) 
in which it defines a mapping between the two set of basis vectors, {|m))}^=i and here. 



III. SQPT IN THE VECTOR DENOTATION 

In the first part of this section, we shall give a brief review of the SQPT theorem which is originally described in 
operator denotation. Then, via the isomorphism developed in the above section, we reexpress the SQPT in the vector 
denotation and show that x matrix can be uniquely decided by the A matrix. 

Suppose that {Djuj^^i is a set of fixed basis operators which satisfy the orthogonal condition Tr[Z)Jj£)y] = 5fiv 
The operators i?™ can be decomposed as E'^ = ^^-^ij- with e™ = Ti[Dj^E"^] the expanding coefficients. Now, 
equation (7) can be rewritten as 

s{p) = J2 D^pDixf,;.. (34) 

Where the coefficients Xfj^-^j^^ 

X,;. = Y,nDiE^]mDtE"^]r, (35) 

m 

are the entries of a rf^ (P matrix x which is Hermitian by definition. 

The central idea of SQPT is to prepare (P linearly independent inputs {Pv}t=i ^^'^ then measure the output e(Pi^). 
In addition, every e(Pi/) can be expressed in terms of a linear combination of basis states. 



£(P,) = ^A^;,P^. 



(36) 
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Since e(P,y) is known from the state tomography, X^^-i, can be determined by standard linear algebraic algorithms. To 
proceed, one may write 

D^PsDi = J2'^^'l^sPf (37) 

7=1 

where are complex numbers which can be determined by standard algorithms from linear algebra given the 
and Pry operators. Combing the last two equations, we have 

E ^jIsX,;. = (38) 

One may take x and A as vectors and ^ as a x rf* matrix with columns indexed by /z; v and rows by 7; i5, and 
rewrite (38) as a set of linear equations 

/3X = A. (39) 

The supper operator x can be determined by inversion of this equation. 

Until now, the general solution of x has not been given. And it has been usually guessed that the solution is 
not unique. In the following argument, we shall show that the vector dentation, which has been developed in above 
section, offers a convenient tool in getting the solution of (39). 

Theorem 8. With x'^ in (8) and the unitary U in (32), a super operator x is defined as 

X = U^X'U, (40) 

its matrix elements, {{iJ^lxW)} , equal with the coefficients Xn;!^ defined in (35). 

Proof. With Theorem 2, we may rewrite (35) as Xti-u = ^rn((^tJ^\^"^))ii^rn,\Dv))- Prom the known definitions, 
X' = Era \E"')){{E'^\ and \D^)) = U\fi)), wc shall get x^;. = ((a^IxI^^)). 

In the vector denotation, the equation (36) can be rewritten \e{P,y))) = ^^.=1 ^i.i;iy\PtJ.))- With the known property 
of the reciprocal states in (24), there should be A^jj, = {{Rf_i\s{P^)) . As it is shown in (15), \s{Pi,)) = X''\P,y)). With 
the aid of theorem 7, we have {{R^l = ((/i|P~^ and \Pv)) = P\i^))- Jointing these results together, we find (36) has 
an equivalent from 

A^;. = {{i,\p-'X^P\u)). (41) 

This result can expressed as the following theorem. 

Theorem 9. With the two super operators , A'^ in (14) and P in (26), a super operator A is defined as 

A = P-^X'P. (42) 

The coefficients Xfj^-^v in (36) should be the matrix elements of it, X^j^-^v = ((/x|A|z^)). 

Note that the two super operators, x and A, are (P x (f matrices in H^^. Via the isomorphism defined in theorem 
3, we could introduce their corresponding vector, |x) and |A) respectively, in the enlarged Hilbert space il®*. Prom 
(40), we could get the relation, x° = UxU\ and express it as 

\x-) = U(^U*\x). (43) 

In the similar way, we may rewrite (42) as 

|A) =P-10P'^|A^) (44) 
Keeping in mind that |A'^) is related to jx'^) with the simple formula in (23), jA'^' — I3'^\x^), we shall arrive at 

\X)=p-^®P'^l3''U®U*\x). (45) 

Let's organize this result as the following theorem. 

Theorem 10. With the given super operators, U , P, and /3'^, a super operator /3 in iJ®* can be defined as 

j3 = P-'^ (S) P'^ P'^U (S) U* , (46) 
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it relates |x) to |A) in the way hke 

|A)=/3|X)- (47) 

The coefBcients, in (38), should be the matrix elements of /3, ^^.'^ = (7; v). 
Proof. At first, we could rewrite (37) in the vector denotation, 

\D,PsDl)) = Y,l3i;:s\P,))- 

7=1 

With the definition of the reciprocal states in (24), there should be 

/3f;-l^,={{R,\D,^ d:\Ps)). 

We could show that the matrix elements of /3 in (46) can be written in the way like 

(7;.5|/3|M;t^) = (|7))((^||P-'®P^/3^C/0C/*||M))((H) 

= {{P-yh)){{6\P^\mn)){{,y\U^) 

= i\R,)){{Psm\D,)){{D^\) 

= {\R^)){{Ps\\D^(g>D:) 

= Tr[{\R,)){{Ps\)^D,<S^D:] 
= {{R^\D^(g> d:\Ps)). 

It should be noted that i^) is equivalent to the vector \ij;kl) with the constraints, as we have assumed in (25), 
/i = (d — l)i + j and v = {d— l)k + 1. Therefore, = \ij)){{kl\ = Cij-ki- In Theorem 3, it has been shown 

Cij:ki should be related to the basic vector \ij;kl). Formally, there is ||/i))((z/|) = |/i;z^). Besides this result, theorem 
4, theorem 7, and (33) have also been used in the derivation above. 

The way of introducing (47) offers a well understanding about the original formula in (39). It clearly shows why 
one may take x a-nd A as vectors and /3 as a rf"* x matrix. At the end of this section, we shall prove that x is 
uniquely determined if A has been given. 

Theorem 11. For the equation (47), |x) is uniquely determined by |A), 

Ix) =r'|A),/3-i =?7toC7T/3^P(^(p-i)T. (48) 

Proof. The equation (47) is nothing else but a system of linear equations. According to the famous Cramer's rule 
[37]: If Ax = b is a system of n linear equations in n unknowns such that det{A) ^ 0, then the system has a unique 
solution. So, if dct(/3) 7^ 0, we may conclude the solution of (47) is unique. To derive the quantity det(/3), we shall use 
the following known results [37]; A. // F and, A are two d^ x d^ matrices, there should be det(FA) = det(F) det(A); 
B. The determinant of a d x d matrix A equals the one of its transpose A"^ , det(j4) = det(j4"'"); and C: If A and B 
are two d x d matrices, there should be det{A B) = det(^) det(i3). 

Noting that both and U are unitary matrices, | det(/3'^)| = | det(C/)| = | det(i7''')| = 1, we shall have 

|det(/3)| = |det(P) dct(p-^)|. 

Considering the fact that we have required the set of operators {Pp}J|^Li to be linearly independent, both the determi- 
nants, det(P) and det(P~^), are nonzero. Therefore, we conclude that the solution of (47) is unique since det(/3) 7^ 0. 
Finally, it should be mentioned that the isomorphism in Theorem 3 is a one-to-one mapping between the matrix and 
its corresponding vector. So, it can be concluded that x is uniquely determined by A. 



IV. THE OPTIMAL SQPT 

Since that equation (48) can be divided into the following three equations 



|x) = f/^0f/^|x'^),|x'^)=^lA^), (49) 
|A) = P-I P'^IA'^), (50) 
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where x A matrices should be given is A*^ has been decided. In the following argument, we shall limit our attention 
to the task of deciding A'^ in experiment. 

As we have shown, the central ideal of SQPT is to prepare (P linearly independent inputs {f^}^=i and measure 

its output e{Pf^). Here, we shall focus on the situation where {f;i}^Li formulate a rank-one positive-operator-valued 
measure (POVM), 

= w^\^^){<i>^\,w^ = Tr[P^], = U (51) 

while our discussion about it, which has presented in the above argument, still holds here. The reason for this 
assumption is clear: The so defined may be easily prepared and its output s{P^) can be directly decided by the 
various protocols of quantum states tomography (QST). It should be noted that e(P^) has always been supposed to 
be known when our general solution in (48) is constructed. In other words, the SQPT scheme itself does not offer a 
way of performing QST. Here, we suppose that the outputs, which lie in the d-dimensional Hilbert space, are decided 
by the linear state tomography [26]. 

To describe the linear state tomography, we should introduce another set of rank-one POVM {Pn}fi^i, 

cP 

P^ = w^\^^){^^\,w^ = Tt[P^], ^ = I,. (52) 
In special, we suppose that the operators in {P^j^Li are linearly independent. Introducing the super operator P, 

p=Y.\p,))m,p\tA) = \p,)), (53) 

The linearly independent condition holds iff dct(P) ^ 0. We use |-R^))) to denote the reciprocal states, 

\R,))^{P-'m), (54) 

which should have the property that 

Y,\R^)){{P^\=Ia^. (55) 
This property directly leads to the linear state reconstruction formula [26] 

HP.))) = Y.\R,)){{PMP'^))) (56) 

or the equivalent form 

HP^))) = = ((P^|e(P.))). (57) 

Noting our task is to decide A'^ instead of a single output state. It can be realized by preparing a series of linear 
independent inputs P. and measuring each output e{Pv) with the projective operator P^. There are altogether 
such coefficients which can be collected as a x 6? data matrix u, 

^=Y. (58) 

With our definitions of the super operators, P in (27) and P in (53), according to the relation (15), we could rewrite 
the coefficient in (57) as 

u;^;. = ((m|(P)U=P|i.)). (59) 
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Put it back into (58), there exists a simple formula connecting oj to A'^ 

LJ = P^XT (60) 

which has an equivalent form 

= {P-yuP-' = ^ u;,AP-y\f^)){W\P-'- (61) 

Recalling our definitions of the reciprocal states, (28) and (54), we shall get one of main results of present work, 

= ^ iv,.,.\R^)){{P^\, (62) 

where the reciprocals states, |-R,j,)) and \Bv)), arc known if both the IC-POVMs, {P^} and {-Pi/}, have been designed 
at the beginning. In the vector denotation, the above equation takes the form 

|A'^)= E ^,;M,)){{PA)- (63) 
Applying the unitary transformation which is defined in (21-22), on both sides of the above equation, we shall get 

|X^)=^=|A^)= 0J^;^\R^®K)- (64) 

In the operator denotation, it should be 

X''= J2 ""i^'-^i^ ® (65) 

fl.U—l 

Before showing which choices of the IC-POVMs are optimal in getting x'^, we shall at first introduce some denotations 
to simplify our discussion. Formally, we define 

= p^(^p;, (66) 

where the subindex x takes value according to the rule 

x={(f - + V. (67) 
One may easily check that{na;}^li is a rank-one IC-POVM in iJ®^, 



l^n, = 1^2, (68) 



and the operators in this set are linearly independent. At the same time, we define the operator Q^-, 

Q^ = R^®Rt. (69) 

Note {Qx\Ily) = Tr[Ql.Uy], it can be easily verified that IQ^) should be the reciprocal states corresponding to Ilia,). 
Therefore, there should be 

d^ 

{QM = |gx)(n,| = ld4. (70) 

x=l 

With theses denotations in hand, one may simplify (64) into a more compact form, 

d* 

|x^) = Eu;(x)|g.) (71) 
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where w{x) = lj^-v holds under the constraint in (67). Timing (11^1 on both sides of the above relation, one may get 

cj{y) = (n,|x^). (72) 

Noting Tr[x'=] = d, we may introduce p^, which is defined in (11), for its normalized state. We rewrite the above two 
equations into such a form 

1 

\pe) = '^P{x)\Qx),P{x) = (n^\pe) = (73) 



x=l 



It's a standard form of linear state reconstruction formula [27]: The unknown state is measured with the rank-one 
IC-POVM {n^l^^i, which Ues in in the Hilbert space, and reconstructed with the canonical dual frame {Qx}t=i- 
From it, we are able to show which choice of 11^: should be optimal by following the general discussion of Scott in 
[27,28]. Suppose that yi, j/2, yjv are the outcomes of measurements on N identical copies of the state pe- The 
estimate for the outcome probabilities is 



1 ^ 

Pe{x) = p{x;yi,y2,...,yN) ■= ■j;;j^S{x,yk), (74) 

k=l 

which gives a linear tomography estimate of ps which obeys the expectation E\p{x)] = p{x). An elementary calcula- 
tions shows hat the expected covariance for A'' samples is 

E[{p{x) - p{x){p{x) -p{x))] = j^{p{x)S{x,y) -p{x)p{y)). 

Now suppose that p(x) are outcomes probabilities for the IC-POVM {11^:}^!^^ in (68). That is , p{x) = {U.x\Pe)- The 
error in our estimate of p^, 



Pe = Pe{yi,y2,-,yN) = ^p{x;yi,...,yN)Q{x). 



x=l 

as measured by the squared Hilbert-Schmidt distance, is 



\pe - Pel? = X] (p(a;) -p(a;))(p(2;) -p(a;))(Qx|<3i/), 

x,j/=l 



which has the expectation 

1 



E[\\p,-p,f] = j^{ip{x)5{x,y)-p{x)p{y))){Q,\Qy) 



= ^{(EpW(Qx|Qx)-Tr[p,2]} 



= ^(A,(g)-Tr[p2]) 

where the quantity Ap{Q) is defined as 

Ap((3) = Ep(a;)(Qx|Qx) 

x=l 

Since that p^ is unknown, it is the quantity Ap{Q) is now of interest. In general, it depends on the detail of p^. This 
dependence can be moved away if we set pe = Peio; u) = uau^Tr[a] =1) and take (Haar) average over all u G u{(P), 



A-s= / dMH(u)VTr[ntW](Q,|Q,) 
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where /Uh is the unit Haar measure. Using Shur's lemma, 

/ did}i{u)uau'^ = -k^lcP- 
Ju{cP) " 

we shall get the averaged quantity 

1 

^r = ^ETr[nx](Qx|Qx). (75) 

Let's express this formula in terms of the IC-POVMs, {P^} and {Pv}, and their corresponding reciprocal states. For 
the IC-POVM {-P/i}jJ=i, we define a super operator 



\p,)){p. 



.^1 • 

Its inverse should be 

d= _ _ 

According to the lemma 17 in [26], there should be 

Tr[(£)-i] > d{d{d + 1) - 1) (76) 

with equality iff {Pfi}'f^=i is a tight rank-one IC-POVM. The same argument is also suitable for the IC-POVM {-P^j^Li 
( iff {-P^I^Li is an IC-POVM ). Introducing the super operator 



with its inverse to be 



^1 Tr[P;] 



(/c*)-i = ^Tr[p;]ii?;))((i?;i. 

Via the similar argument above, there should be 

Tr[(/C*)-^] > d(d(d + l) - 1). 
For the result in (75), one may easily find that 

A-s = lTr[(£)-i]Tr[(/C*)-i] 

and the optimal value of it should be 

A7s = (d(d+1)-1)2 (77) 

which can be achieved at if the both the IC-POVMs, {P/i}^Li and {-P/JIJJLd have been chosen to be the tight rank-one 
IC-POVM. 

Let's organize the content of the present section in following way: To perform the SQPT for a d— dimensional 
system, we prepare a set of d? linearly independent inputs {Pi,}f,^i, which are supposed to formulate an IC-POVM, 
and measure each output e{Pv) with anther set of IC-POVM {P^}jJ=i- All the experiment data arc presented 
with the matrix uj defined in (58). Within our present scheme, in which the SQPT theory is represented in the 
vector denotation, we could give the way of reconstructing A'^ in (62) and the one for x'^ in (64). As an interesting 
observation, wc find that equation (64), where \x'^) can be viewed as an unknown states in iJ®^, is just a special 
case of the linear state reconstruction formula discussed by Scott in [27]. This observation makes it possible for us 
to get the optimal IC-POVMs, {Pv}'l^i and {P^j^^i, by following the general linear state tomography theory. The 
conclusion, in which {P^j^L^ should be a tight rank-one IC-POVM, can be viewed as an known result in [26]. In 

addition to it, we show that the conjugate of inputs {Pu}'l^x should also be prepared as a tight rank-one IC-POVM. 
This two requirements should make the statical error in reconstructing x'^ to be minimal. 
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V. THE SYMMETRIC INFORMATIONALLY COMPLETE POVM 

Under the condition of linearly independence, the tight rank-one IC-POVM is unique [26]: It's the so called 
symmetric informationally complete POVM (SIC-POVM) introduced in [28] 

1 1 

P, = ^\%){^,\^^%,Y.P,=U (78) 
Where the normahzed states |$^) have the following property [28] 

K*.i*.)i^ = ^rrf- ^^^^ 

From this definition, we see that {-P^j^Li should also be a SIC-POVM. In the following argument, we will suppose 
that input and measurement IC-POVMs are the same one, 

P = P 

The SIC-POVM defined above should also be the 2-design [26] 

d 

T.P^'^P^ = TTTTT^ E ® + \3){k\ ® \k){j\. 

Let's take it as our starting point here to derive the reciprocal (canonical dual frame) operators needed in our 
reconstructing formula (65). Noting is a projective operator which has the property that P* = P^ . Therefore, 
performing the partial transposition on above equation, we shall get 



d'^ ^ d 



(80) 



In the vector denotation, it can be rewritten as 

T.\P^^P*^)=d{lTT)^ lifc;ifc) + |ii;fcfc) 

p=l j,k=l 

in which wc have use our denotation in (17), kl) = (S) 1^0)- Performing the unitary transformation, 13" in (20), 
on both sides of it, we shall get 

^ d 

Noting in our definition of the isomorphism in theorem 3, the basis vector kl) has a one-to-one mapping with the 
operator Cij-^ki in (12). Recalling our definition of the supper operator in (30), we shall get the formula 

Using the relations that Ej,fe=i \k){k\ = 1^2 and J2'j,k=i b')(^l <^ = we shall get 

^=^(I.^ + |I.))((I.|) (81) 

with its inverse to be 

:F-^=d{d+l)Id^-d\ld)){{Id\. (82) 
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With equation (31), we find the reciprocal states 

\R^))=J^-'\P^))=d{d+l)\P^))-dTr[PM 
Prom our definition of in (65), we shall get a compact form of i?^, 

i?^ = (rf+l)$^-ld, (83) 

which can be used to reconstruct x'^ in (65). It's meaningful to calculate another super operator 



(84) 



with it's inverse to be 

/C-i = (d+l)Irf2-|Id))((Id|. (85) 

Formally, it can be rewritten as the mixture of a set of basis vectors {|eju))}^Li, which satisfy the constraints that 
= 7^ Ej=i lii)) and ((e^|e^)) = 5^^, 

lC-' = \e,)){{e,\ + {d+l)J2K)){{e,\. (86) 

Obviously, its eigenvalues should Ai = 1 and = d + 1 for /i 7^ 1 . The trace of it should be 

Tr[/C-i] = l + (d+ l){d^ - 1) = d{d{d + 1) - 1) (87) 

in which the lower bound of (76) is achieved here. By the way, one may check that /C* = /C, the above bound also 
holds for the SIC-POVM {P*}. 

With the assumption that both the inputs and measurement operators are the SIC-POVM defined in (78), we may 
rewrite the coefficients uijj,-^ in (57) as 

'^m;- = ^TV[$t£($.)]. (88) 
Putting (83) back into (65), we shall get our reconstruction formula for x'^, 

= -Id- + "^t^-Aid + 1)'*^ <»K-{d+ 1)*^ Id], (89) 
1 

in which we have used the relations ^t^^v = d and Y^^^i = 1/c?. If the quantum channel is unital, where 

e{ld) = Id, the above formula can be simplified as 

d= 

X" = -{d + 2)1^2 + {d+ if J2 ® (90) 

In the end of this section, we shall offer an example to show how our reconstruction formula works. Let's consider 
the depolarizing channel, where an arbitrary state $ should have its output to be 

£($) = (l-g)$ + ^Irf, (91) 

with q a positive parameter taking values form to 1. One may check it's an unital quantum channel since £(Id) = Id- 
Via a simple calculation, we find 

Tr[$t £($,)] = -^-T^ V + 



d+1 d{d+l)' 
Jointing it with equations (80) and (90), we find that x'^ is now a Werner- type state. 



X" = dpwevnev, PWeruev = (1 - g)|S'+))((S'+| + ^1^2, (92) 



where \S^)) is the maximally entangled state defined in Theorem 1. 
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VI. DISCUSSION 

In the ancilla-assisted QPT, one may prepare a maximally entangled state \S+)) = \ ^ input 

with its corresponding output defined as 

p, = £0ld(|5+))((5+|). 

If this output has been decided with the quantum state tomography, the information of the quantum process, which 
is now described by the the matrix x"^ = X]m should also be known as x'^ = dpe- 

As we have shown, our present protocol is equivalent to deciding pe with a product SIC-POVM. In [27], Scott has 
demonstrated that the so-called unitary 2- designs are the optimal choice for the ancilla-assisted QPT in getting p^. 
For the expectation 

^[||Pe-/5e|P] = ^(A,(Q)-Tr[p,n, 
the unitary 2-designs should offer a optimal value 

A°P*(Q)=d4-d2 + J_, (93) 

for the general quantum channel ( an equivalent result can be also found in [38]). There is a simple inequality between 
this optimal value and our result A^^^ in (77): 

When d, the dimension of the system, is increasing, the ratio above should approach unit. A related discussion about 
the product SIC-POVM can be also found in [39], where the authors shows that there is only a marginal efficiency 
advantage of the joint SIC-POVM over the product IC-POVM in the state tomography. 

In the present work, we have assumed that the outputs are measured with a set of linear independent operators. This 
assumption make it possible for us to introduce the Chefles' reciprocal states which is convenient for our calculation. 
However, a set of nonlinear independent IC-POVMs can be also used in the quantum state tomography. How to 
generalize our scheme for such situations, should be studied in future. 

Finally, let's end our work with a short conclusion. At first, we developed a self-consistent scheme, in which 
bounded linear operators are presented by vectors, for constructing the set of linear equations in order to get the 
matrix containing the complete information about the qiiantiim process. We proved that the solution of this set 
of linear equations is unique. In second, letting the inputs and the measurements be prepared by two sets of the 
rank-one positive-operator- valued measures [POVMs] , where each POVM is supposed to be linearly independent and 
informationally complete (IC), we observe that SQPT now is equivalent to deciding a unknown state with a set of 
product IC-POVM in the d^-dimcnsional Hilbert space. Following the general linear state tomography theory, we 
show that the product symmetric IC-POVM should minimize the mean-square Hilbert-Schmidt distance between 
the estimator and the true states. So, an optimal SQPT can be realized by preparing both the inputs and the 
measurements as the symmetric IC-POVM. 
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